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ABSTRACT 

We report new observations of the galaxy UGC8802 obtained through GASS, the GALEX 
Arecibo SDSS Survey, which show this galaxy to be in a remarkable evolutionary state. UGC8802 
(GASS35981) is a disk galaxy with stellar mass M* =2x 10 10 M Q which appears to contain an addi- 
tional 2.1 x 10 10 M Q of HI gas. New millimeter observations with the IRAM 30m telescope indicate 
a molecular gas mass only a tenth this large. Using deep long-slit spectroscopy, we examine the spa- 
tially resolved star formation rate and metallicity profiles of GASS35981 for clues to its history. We 
find that the star formation surface density in this galaxy is low (Esfr. = 0.003M Q yr~ x kpc~ 2 ) and 
that the star formation is spread remarkably evenly across the galaxy. The low molecular gas masses 
measured in our three IRAM pointings are largely consistent with the total star formation measured 
within the same apertures. Our MMT long-slit spectrum reveals a sharp drop in metallicity in the 
outer disk of GASS35981. The ratio of current star formation rate to existing stellar mass surface 
density in the outer disk is extremely high, implying that all the stars must have formed within the 
past ~1 Gyr. At current star formation rates, however, GASS35981 will not consume its HI reservoir 
for another 5-7 Gyr. Despite its exceptionally large gas fraction for a galaxy this massive, GASS35981 
has a regular rotation curve and exhibits no sign of a recent interaction or merger. We speculate 
that GASS35981 may have acquired its gas directly from the inter-galactic medium, and that it and 
other similar galaxies identified in the GASS survey may provide rare local glimpses of gas accretion 
processes that were more common during the prime epoch of disk galaxy formation at z ~ 1. 
Subject headings: galaxies: individual (UGC8802) - galaxies: star formation - galaxies: evolution - 
galaxies: ISM - galaxies: stellar content - galaxies: kinematics and dynamics 



1. INTRODUCTION 

One of the largest gaps in our understanding of how 
galaxies form and evolve is the question of how gas — 
the raw material for star formation — flows into and out 
of galaxies, and how these flows regulate star formation 
in these systems. Although ACDM simulations make 
specific predictions for how structure in the dark matter 
assembles through hierarchical clustering, the assembly 
of the visible, baryonic components of galaxies is still a 
subject of considerable controversy. 

The currently favored theoretical picture is that cold 
gas flows along filaments into the centers of assembling 
dark matter halos at high redshifts (z > 2), and this 
process can build massive galaxies rapidly at early times 
(e.g., Genel et al 2008; Dekel et al 2009). At lower 
redshift, massive spiral galaxies are thought to accrete 
gas much more slowly, at a rate of less than a few solar 
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masses per year, from a surrounding hot corona (Mailer 
& Bullock 2004; Kauffmann et al 2006; Peek, Putman & 
Sommer-Larsen 2008; Binney, Nipoti & Fraternali 2009). 
Gas may also continue to accrete through major or minor 
mergers with other galaxies, although this is not believed 
to be the primary way in which ongoing star formation 
in local spiral galaxies is fuelled at present (Sancisi et 
al. 2008). Outflows, whether from star-formation or 
AGN-related feedback mechanisms, are another piece of 
the puzzle, and their role in shutting down star forma- 
tion, limiting the efficiency with which stars are able to 
form, or even providing fuel for ongoing star formation 
in the form of recycled wind material (Oppenheimer et 
al. 2009), is similarly debated. 

In order to learn more about cold gas in nearby galax- 
ies, we are carrying out the GALEX Arecibo SDSS Sur- 
vey (GASSJ3 (Catinella et al. 2010). GASS is designed 
to measure the neutral hydrogen content of a represen- 
tative sample of ~ 1000 galaxies uniformly selected from 
the SDSS spectroscopic and GALEX imaging surveys, 
with stellar masses in the range 10 10 — 1O 115 M0 and 
redshifts in the range 0.025 < z < 0.05. As GASS obser- 
vations are designed to detect HI down to a gas-fraction 
limit of 1.5%, the full GASS sample will be the first HI 
survey able to place meaningful, unbiased constraints on 
the atomic gas reservoirs that may contribute to future 
growth in massive galaxies. 

We are also pursuing a companion project on the 

8 http://www.mpa-garching.mpg.de/GASS 



IRAM 30m telescope, COLD GASlfi which will obtain 
accurate and homogeneous molecular gas masses for a 
subset of ~ 300 galaxies from the GASS sample. These 
data will allow us to characterize the balance between 
atomic and molecular gas in the galaxies in our sample, 
and understand the physical processes that determine 
how the condensed baryons are partitioned into stars, 
HI and H2 in the local Universe. 

In this paper, we report on UGC8802, an extraordi- 
nary galaxy blindly selected for inclusion in the GASS 
sample (under the catalog name GASS35981, used here- 
after), which contains a reservoir of HI > 10 10 M©, at 
least equal in mass to the galaxy's entire stellar con- 
tent. This galaxy contains less than one tenth this mass 
in H2 and also has a rather modest star formation rate 
(SFR). We describe the spectroscopic follow-up that en- 
ables us to conclude that the outer disk of this galaxy is 
currently forming from gas that has likely accreted from 
the external environment. In the following, we adopt a 
standard ACDM cosmology with Ho = 70km s _1 Mpc -1 , 
_. m = 0.3 and ^a = 0.7. 

2. OBSERVATIONS 

GASS35981 (also SDSS J135308.36+354250.5, in addi- 
tion to UGC8802), was selected for inclusion in the GASS 
parent sample because it has photometry from both 
SDSS and GALEX, is located in a region of sky accessible 
to Arecibo, and has stellar mass of M* = 2 x 10 10 M© and 
redshift of z = 0.0411 that fit into our targeted range. 
GASS35981 has pre-existing HI observations available 
from the Cornell HI Digital Archive (Springob et al. 
2005). The mass of HI in GASS35981 is estimated from 
the line flux to be 2.1 x 1O 1O M0, and the rest-frame ve- 
locity width of the HI line is W 50 = 360±25 km s" 1 (Fig- 
ure [TJ. We display the HI archive spectrum in Figure Q] 
We note that there is no evidence for contamination 
from possible companion galaxies within the 3.5' Arecibo 
beam: one nearby galaxy has a spectroscopic redshift of 
0.14, and two additional faint companions have SDSS 
photometric redshifts consistent with z — 0.14. 

The lower panel of Figure ED shows that GASS35981 
(blue star) lies near the extreme end of HI fractions ob- 
served by GASS. Indeed, its gas fraction is comparable 
to the highest values measured for all galaxies in this 
stellar mass range (e.g., Giovanelli et al. 2007). Indeed, 
even in Hi-selected samples, which are biased towards ob- 
jects like GASS35981, galaxies in this stellar mass range 
(M* > 10 10 M Q ) with such high gas fractions are quite 
uncommon (e.g., the "HI Giants" described by Garcia- 
Appadoo et al. 2009). 

It is also clear that GASS35981 lies significantly above 
the best-fit 'gas fundamental plane' (dotted line) relat- 
ing stellar mass surface density, NUV-r color, and gas 
fraction, as described in Catinella et al. (2010). Because 
this galaxy was an interesting outlier and a presumed 
easy target, GASS35981 was selected for inclusion in our 
initial COLD GASS pilot program to obtain molecular 
gas measurements with the IRAM 30m telescope. 

Observations of GASS35981 in the J=l-0 rotational 
transition of CO were made at 3mm with the IRAM 30m 
telescope, in three different pointings: one at the galaxy 
center, and one each to the north and south, one bcam- 

9 http://www.mpa-garching.mpg.de/COLD_GASS 
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Figure 1. Top: HI spectrum of GASS35981 from the Cornell HI 
Digital Archive (Springob et al. 2005). The solid orange line in- 
dicates the measured Hi line center, and the dashed lines mark 
the velocity width of the spectrum. Bottom: GASS35981 on the 
'gas fundamental plane' presented in Catinella et al. (2010), re- 
lating HI mass fractions to stellar mass surface density (/J.*) and 
NUV-r color. The best-fit values of the coefficients a, b, and c 
used in constructing the plane are: a = —0.332, b = —0.240, and 
c = 2.856. GASS35981 is marked with a blue star. Red circles in- 
dicate Hi-detected galaxies from GASS, and green triangles mark 
upper limits on gas fraction. Blue circles denote a set of galax- 
ies randomly drawn from previous HI surveys, and added to the 
GASS galaxies in the proper proportion to generate a statistically 
unbiased sample (see Catinella et al. (2010)). 

width (22") away along the galaxy major axis. Data 
were taken in June and August 2009, using the WILMA 
and 4MHz backends simultaneously to record the data, 
and the CLASS_3 software to process them. Individual 
scans were examined, and a linear baseline subtracted 
from each of them. After rejection of scans with unstable 
baselines due to, e.g., poor atmospheric conditions, the 
data were combined and binned to a spectral resolution 
of 21 kms" 1 . 

The CO line is detected in the central pointing with 
S/N= 5.8, for an integrated line flux of T* = 0.43 ± 
0.07K km s _1 within a 400 km s _1 -wide window. We set 
upper limits at 0.29 and 0.43 K km s _1 for the North 
and South offset pointing, respectively. Adopting a con- 
version factor of Xco = 4.4M©/L' (where L' has units of 

10 http://www.iram.fr/IRAMFR/GILDAS 
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Figure 2. CO spectra for all three IRAM pointings, as labeled on the figure. The red lines mark the optical redshift of GASS35981, the 
blue dashed lines show the velocity centroid of the central CO emission, and the dotted lines represent the width of the window function 
that was used to measure the flux in the CO lines, covering an area marked by grey shading under the spectrum. In the central pointing, 
the peak at ~ 300 km s — 1 is believed to be noise, as it does not appear in both polarizations; it therefore is excluded from our measurement. 



K km s 1 pc 2 ), these fluxes correspond to an H2 mass 
of 8.8 x 10 8 M Q in the central 22" pointing, and upper 
limits of 5.9 and 8.8 x 10 8 M Q in the offset pointings. Us- 
ing an aperture correction for the 22" beam of the IRAM 
telescope based on resolved CO maps of nearby galaxies, 
we estimate a total H2 mass of 1.45 x 10 9 M for the 
galaxy, based on the detection in the central pointing. 
Following the prescription of Springob et al. (2005), we 
measure a line width of W50= 335 ± 20km s^ 1 . This 
corresponds to a rest- frame width of 321 km s _1 . The 
spectra of all three pointings are shown in Figure [2] 

Follow-up long-slit spectroscopy of GASS35981 was ob- 
tained on 20 November, 2009, using the Blue Channel 
Spectrograph on the 6.5m MMT telescope on Mt. Hop- 
kins, AZ. The spectrum was obtained with a slit of width 
1.25"oriented to PA= 13.4° on the sky, such that the slit 
runs along the major axis of the galaxy, as indicated 
in the second panel of Figure |3] GASS35981 was ob- 
served in 2x600s exposures with the 500 line grating. The 
spectrum covers a wavelength range of ~ 3900 — 7000A 
at a spectral resolution of ~ 4 A FWHM, equivalent to 
a ~ 75km s _1 in the rest frame of the galaxy. The slit 
samples the galaxy spatially with 0"3 pixels (equivalent 
to 0.25 kpc in physical units), over a total slit length of 
150", which is much larger than the galaxy itself. 

Data were reduced in IDL with routines from the pub- 
licly available Low-Redux packageFI maintained by X. 
Prochaska, which itself is an adaptation of the Princeton 
SDSS reduction codeP^I to more general long-slit reduc- 
tions. The code performs standard biasing, flat-fielding, 
cosmic-ray rejection, and sky estimation on each ex- 
posure. We then co-add the sky-subtracted exposures 
through a custom-written routine that verifies and ad- 
justs the alignment of exposures before co-addition. Flux 
calibration was achieved via observation of the spec- 
trophotometric standard BD+17 4708. We note that 
we have not attempted to apply any correction for at- 
mospheric dispersion, even though GASS35981 was ob- 
served at an airmass of 1.8. The expected atmospheric 
dispersion from blue to red end for our configuration, as 

11 http://www.ucolick.org/~xavier/LowRedux/ 

12 http://spectro.princeton.edu 



estimated from Filippenko & Greenstein (1984), is com- 
parable to the size of our slit. More importantly, it is 
much smaller than the typical spatial regions we aver- 
age over in the analysis presented below. A cut-out of 
the fully-reduced spatially-resolved spectrum centered on 
the Ha emission line, including only that portion of the 
slit where significant flux from the galaxy is detected, is 
displayed in Figure |3l 

3. ANALYSIS 

To understand why GASS35981 contains such a high 
abundance of HI, but little molecular gas, we analyze our 
long-slit spectrum with an eye toward extracting infor- 
mation on the dynamics of the galaxy and the properties 
of its stellar population — including current star forma- 
tion rate, stellar population age and metallicity, and dust 
properties — as a function of radius across the galaxy. We 
likewise re-analyze the SDSS ugriz and GALEX FUV and 
NUV images of GASS35981 to obtain independent esti- 
mates of the radial gradients in colors and star formation 
histories, as described below. 

3.1. Rotation Curve 

To measure the rotation curve of GASS35981, we first 
bin our spectrum spatially to achieve a signal to noise 
of at least 6.0 (A -1 ) in each bin, the minimum needed 
for a reliable velocity estimate. We define spatial bins by 
working outward from the galaxy center (where S/N is 
high even without binning), co-adding rows one at a time 
until the minimum S/N is reached, or until adding rows 
fails to increase the S/N. Since strong Ha emission is de- 
tected in GASS35981 out to much larger radius than the 
continuum flux (see Figure [3|), we also include several 
additional bins that correspond to knots of bright Ha 
emission in the spectrum. For each binned spectrum, 
we determine the effective spatial position by calculat- 
ing the luminosity- weighted average radius of all spatial 
positions that entered the co-add. 

We then measure the radial velocity as a function of 
radius in two ways: 1) by cross-correlating each spec- 
trum against galaxy templates to determine the velocity, 
2) by fitting directly for the centroids of the Ha emission 
line. The method utilizing Ha allows us to extend our 
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Figure 3. Top: GALEX UV image of GASS35981 with FUV 
in blue and NUV in green, smoothed by a Gaussian with FWHM 
3.5", slightly smaller than the GALEX psf of 4.5". The dimensions 
of the image are 87" wide by 42" tall. Upper Middle: Color SDSS 
image of GASS35981 constructed from gri images, and displayed at 
the same scale and size as the GALEX image above. Orientation 
and width of the spectroscopic slit are overlayed in red, with north 
indicated by the compass rose. White circles overlayed indicate 
the three IRAM CO pointings, with circle diameter equal to the 
22"FWHM IRAM beam size. Lower Middle: A cutout of the 
spectrum around the region of the Ha emission line is shown at 
the same spatial scale along the x-axis. Bright Ho knots can be 
clearly associated with faint blue regions of the outer galaxy that 
intersect the slit. Vertical red lines mark out the twelve spatial 
regions over which we have co-added our spectrum. All spectro- 
scopic measurements in this paper are measured on these twelve 
spectra. Bottom: NUV-r color profile of GASS35981 extracted 
from a region overlapping our slit. Note that the X-axis arcsccond 
scale applies to all four panels of this figure. 

rotation curves to larger radii. The two methods yield 
consistent curves at radii where we can fit templates and 
measure emission line centroids. The resulting rotation 
curve is plotted in Figure |4l not corrected in any way 
for the inclination of the galaxy. The solid line plot- 
ted over the measured curve is the best-fitting rotation 
curve of the form V(R) = V M AxR/(R a + R") 1/a + AV 
(B0hm et al. 2004; Moran et al. 2007), where R is 
the radius, a and R s are free parameters that govern 
the shape of the rotation curve and its turn-over, and 
AV is the offset of the galaxy's central velocity from the 
redshift obtained from the SDSS spectrum (also left as 
a free parameter). We measure a circular velocity for 
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Figure 4. Rotation curve of GASS35981 measured primarily from 
cross-correlation of the spectrum to templates, supplemented by 
measurements of the Ho centroid for the two outermost bins to 
each side. Solid line indicates the best-fitting model of the form 
described in the text, which is symmetric to positive and negative 
positions. 

GASS35981 of 2F M ax = 318 ± 10 km s" 1 . This is con- 
sistent with the value obtained from the CO spectrum 
(VK = 321±10kms- 1 ). 

Table 1 summarizes some of the basic properties of 
GASS35981, including the masses and line widths of 
the various components. The agreement between CO 
and Ha widths is not surprising because at the ra- 
dius of the IRAM beam (11 arcsec), the optical veloc- 
ity has already reached its maximum value. It is inter- 
esting, however, that the rest-frame HI velocity width, 
W — 360 ± 25 km s _1 , appears to be somewhat larger 
than the estimates from both CO and Ha. Ej Since the 
optical rotation curve of GASS35981 is fiat, one might 
have expected a better agreement between the HI and 
Ha widths. 

Even though the significance of the difference is 
marginal (smaller than 2a), it is interesting to comment 
on possible causes of the effect. Based on a statisti- 
cal analysis of a large data set with both Ha and HI 
spectroscopy, Catinella et al. (2007) find a systematic 
difference between optical and HI velocities of normal 
spiral galaxies with flat rotation curves of a factor 1.06 
(see their figure 6). After applying a 6.5 km s _1 tur- 
bulent motion correction (as in Catinella et al. 2007), 
and accounting for such a scaling factor, the HI width 
of GASS35981 is 333 km s" 1 , which is in slightly bet- 
ter agreement with the optical rotation curve and CO 
linewidth. The residual difference might be explained by 
uncertain statistical and turbulent motion corrections, 
or it might provide evidence for accretion of HI at large 
radii. It is interesting that the HI profile shown in Fig- 

13 We note that the velocity widths are not deprojected to edge- 
on (the galaxy's inclination is 63.25 degrees from face-on, estimated 
as cos(i) = b/a = 0.45). 
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Figure 5. Example of absorption template (top) and emission line (bottom) fits to the observed spectrum of GASS35981, here shown 
for a spectrum extracted 2-3kpc from the center of the galaxy. In the top panel, the spectrum is plotted in black with regions that were 
masked for the fitting indicated in green. The best-fitting model spectrum is over-plotted in red. In the bottom two panels, two cutouts 
of the spectrum near the H/3 and Ha emission lines are replotted with the best-fit model subtracted off. Our Gaussian fits to key emission 
lines are then over-plotted in red. 



ure [T] is noticeably asymmetric, 
this point in the final discussion. 



We will come back to 



3.2. Spatial Binning of the Spectrum 

We now correct each row in the unbinned 2D spectrum 
to a common rest-frame using the rotation curve in Fig- 
ure |U This procedure is necessary if we wish to obtain 
accurate measurements of spectral lines. We then re-bin 
the spectrum to a slightly higher targeted S/N of ~ 8. 
To each side, the two bins at highest radii are chosen by 
hand, because there is no clear detection of the stellar 
continuum this far out. As can be seen in Figure El neb- 
ular emission is still clearly detected in these regions, so 
we choose the binning to correspond to clear transitions 
between the bright or faint knots of Ha. We end up with 
a total of 12 radial bins, and we mark the limits of each 
with red lines on the spectrum in Figure El We use this 
re-binned spectrum in our subsequent analysis. 

3.3. Continuum & Emission Line Fitting 

We employ a modified version of the technique of 
Brinchmann et al (2004) to measure the strengths of 
key emission and absorption lines as a function of ra- 
dius across GASS35981. The spectrum from each of our 
12 spatial positions is first fit to a linear combination 
of templates drawn from Bruzual & Chariot (2003) sin- 
gle stellar population models, masking regions of the 



spectrum occupied by common emission lines. Best- 
fit model spectra are generated individually for both a 
low- and high-metallicity set of Bruzual & Chariot mod- 
els (0.02Z Q and Z Q , respectively), and we retain the 
one with lowest reduced \ 2 °f the fit. Near the center 
of the galaxy, these fits unambiguously favor the solar- 
metallicity model. Low-metallicity templates provide a 
slightly better fit in each of the two outermost bins where 
we can measure the stellar continuum ( R ~ 15 kpc on 
each side). However, we caution that the S/N in the 
continuum is < 5 here, and so the best-fitting metallicity 
is unlikely to be a reliable indicator. The results of an 
example fit are shown in Figure [5] 

Next, we subtract the best- fitting stellar contin- 
uum model from the measured spectrum, creating an 
emission-line only spectrum where the Balmer emission 
lines can be measured free of contamination from the un- 
derlying stellar absorption. In the two outermost bins, 
only nebular emission is detected and no continuum fit- 
ting is possible. In these cases, we fit a low-order polyno- 
mial to the spectrum to correct for small imperfections in 
our sky subtraction, which arise when co-adding across a 
large portion of the slit. We then measure emission lines 
using the polynomial-subtracted spectrum. 

We fit a Gaussian function to the emission lines, with 
the width of the Gaussian constrained to a single value 
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a Springob et al. (2005) 

b Rest-frame line width 

c SDSS r-band isophotal major axis radius (25mag/" ) 

d SDSS r-band Petrosian half-light radius 



for all lines in a given spectrum. The positions of the 
line centroids are constrained to their rest- wavelengths. 
The only free parameters in the fit are the amplitude of 
each line and an overall velocity offset term. For each 
line, we also recalculate the continuum level in a small 
region within 50A of each fitted line, and subtract this 
small residual from the template fits to further refine the 
measured emission line flux. In addition to the Balmer 
lines Ha, H/3, and H7, used to estimate dust extinction 
and star formation rate (see below) , we also measure the 
forbidden lines [OIII] 5007, [Nil] 6548 and 6584, and 
[SII] 6717 and 6731 to measure metallicity across the 
galaxy. 

We limit our analysis of stellar absorption features to 
the .D4000 n index, which measures the strength of the 
4000A break. This index is an indicator of stellar pop- 
ulation age (see, e.g., Kauffmann ct al. 2003). We use 
Z?4000 n values measured from the best-fitting absorption 
line templates, rather than from the spectrum itself. The 
two agree well for spectral bins with high S/N. We choose 
the model measurement as we expect it to yield a more 
reliable estimate in bins where S/N in the continuum is 
only ~ 3 — 5, because the model fit utilizes information 
from the entire spectrum. 

We estimate the dust extinction within the nebular 
gas by calculating the Balmer decrement, which we de- 
fine as the ratio of Ha/H/3 to the case B recombina- 
tion ratio of 2.87 (Osterbrock 1989). We adopt the for- 
mula A v = 1.9655i?ylog(Ha/H/3/2.87), where we as- 
sume Ry = 3.1 and adopt the Calzetti (2001) extinction 
curve. We further refine the estimates of Ay by com- 
paring to the H7 line, which is only strong enough to be 
measured in emission in a few cases. However, even as 
an upper limit, H7 provides a useful constraint. If, after 
correcting all three lines for dust based on our original 
Ay estimate, the ratios H7/HQ; and/or H7/H/3 are incon- 
sistent with the expected values of 0.474 and 0.166, re- 
spectively, we adjust Ay upward or downward — staying 
within the la errors on the original measurement — to im- 
prove the agreement with these values as far as possible. 

We note that all SFRs reported below have been cal- 
culated after correcting fluxes for extinction using the 
Balmer decrement. Likewise, our equivalent widths have 
been corrected for dust using our measured Balmer decre- 
ments and assuming E(B-V) of the stellar light is ~ 0.44 
E(B-V) of the gas, as in Calzetti (2001). 

3.4. Colors and Star Formation Rates from the 
Photometry 



We measure the UV/optical color profiles of 
GASS35981 using SDSS and GALEX photometry in two 
ways. First, we measure fluxes in a 5" wide strip running 
along the major axis of the galaxy. Our chosen aperture 
lies on top of, but is wider than the spectroscopic slit. To 
ensure that the fluxes measured using the SDSS images 
are well-matched to those measured from the GALEX 
images, we convolve the SDSS images with a Gaussian 
of FWHM 4.5"(the size of the GALEX PSF). The re- 
sulting color profile is displayed in the bottom panel of 
Figured 

For more precise comparison to the spectroscopy (but 
at the expense of lower S/N in lower surface brightness 
regions), we also measured colors through apertures with 
widths identical to the slit segments containing each of 
our twelve spectral bins. These apertures are smaller 
than the GALEX PSF, so we scale the SDSS w-band flux 
measured in the smaller apertures by the ratio of FUV 
or NUV to u measured in the wider aperture at the same 
position along the slit. We use these measurements to es- 
timate SFR and stellar mass surface densities by fitting 
the full UV through optical SEDs to population synthe- 
sis models as described in Wang et al. (2010, in prepa- 
ration). Our method is very similar to that described in 
Salimet al. (2007). 

To optimize the fits, we constrain the internal extinc- 
tion of the galaxies in the library of models to be consis- 
tent with the measured Balmer-decrement values for this 
galaxy. Assuming that the attenuation of the starlight 
is 0.44 times the extinction measured in nebular gas 
(Calzetti 2001), we obtain t v — 0.5 averaged over our 
12 spectral bins, with a dispersion of ±0.25. We use this 
as a prior in our parameter estimation. 

4. RESULTS 

In this section, we analyse star formation rates, stellar 
population ages and metallicities of the nebular gas as a 
function of position in GASS35981. We show that star 
formation is spread fairly evenly across the disk of the 
galaxy, but the ratio of current star formation to stel- 
lar mass surface density in the outskirts of the galaxy is 
very high. This suggests the outer disk has been formed 
recently. We then show that the modest H2 content im- 
plied by the CO observations matches expectations from 
the star formation rate profile of the galaxy. We also 
examine the metallicity profiles, finding a sharp drop in 
metallicity that corroborates our conclusion that stars in 
the outer disk were formed recently. 

4.1. Star Formation in GASS35981 

After correcting Ha luminosities for dust, we measure 
star formation rates using the equation in Meurer et 
al. (2009): SFR (M yr- 1 )=L Ha /(6.93 x 10 33 W), cor- 
rected to a Kroupa (2001) IMF. Then, by dividing by 
the area of the galaxy covered by each portion of the 
slit (1.25" x 0.3"N, where N is the number of individual 
rows that went into each co-added spectrum), we esti- 
mate the SFR surface density as a function of position 
across GASS35981. This is shown in Figure EJ left, from 
which it is clear that the SFR surface density is remark- 
ably uniform across the galaxy. This result is confirmed 
by the SFR surface densities estimated from the photom- 
etry alone (§ 3.3), over-plotted as red circles on Figure|6l 
which match the values estimated using the spectroscopy 
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Figure 6. Left: Star formation rate surface density, in Mq yr —1 kpc -2 , as a function of radius for each spatial position across GASS35981. 
Black points indicate SFR surface density measured directly from dust-corrected Ha luminosity, while red circles are estimated from fits 
to the photometry as described in the text. Current star formation appears to be evenly spread across the galaxy, with a scatter of only 
0.2 dex. Right: EW(Ha) as a function of position (black points and lower limits), which can be expressed as the ratio of current to 
past-averaged star formation, b, as described in the text. Red circles indicate b values estimated directly from the photometric fits. 

pends on the past star formation history of the galaxy. 
An empirical relation between the two quantities can be 
derived for galaxies in the SDSS using the star forma- 
tion rates, stellar masses and Ha equivalent widths pub- 
lished in the SDSS MPA/JHU value added catalogO 
The 6-parameter is estimated from the SFRs and stellar 
masses using the formulation described in Brinchmann 
et al. (2004). We use this relation to transform between 
EQW(Ha), which is plotted on the left axis of Figure [HI 
to b, which is plotted on the right axis of the same figure. 

Figure [6] shows that the timescale for building up the 
current stellar mass surface density at the current star 
formation rate decreases strongly as a function of ra- 
dius. At radii R > 10 kpc we find that b > 1, indicating 
that the SFR is currently higher than it has been over 
ic most of the galaxy's history (Rocha-Pinto et al. 2000). 

b reaches values in excess of 10 in the very outer regions 
of the galaxy, indicating that all the stars in the outer 
disk could have been formed over a timescale as small as 
a few hundred Myrs. 

At large radii, the Ha emission line is very strong, but 
there is no detection of the underlying continuum, yield- 
ing only a lower limit on the EW. We also estimated b 
directly from our photometric fits, and these values are 
over-plotted as red circles in Figure |U Overall, the re- 
sults agree well, though there are two points where the 
estimated b is significantly lower than the limits implied 
by the spectroscopy. Reasons for the discrepancy could 
include small errors in sky subtraction on the spectrum 
leading us to underestimate the 3-er limit on the contin- 
uum flux, which would subsequently cause an overesti- 
mate of b. Alternatively, since the discrepant points are 
in regions where b appears to have a strong spatial gra- 
dient, it may be that seeing and/or resolution effects are 
limiting how well we can match spatial extraction regions 
between spectrum and photometry, and the discrepancy 
simply reflects that spatial variation. In either case, both 
the spectrum and photometry agree that b > 1 all across 
the outer disk. 
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Figure 7. D4000 n vs radius, with the index estimated from the 
best-fitting continuum model of the observed spectrum. In a simple 
single stellar population model from Bruzual & Chariot (2003), 
-D4000 n = 1.3 indicates a stellar population age of 500 Myr, while 
-D4000 n = 1.7 occurs at an age of ~ 3 Gyr. 

very well. The star formation rate surface density profile 
of GASS35981 is quite different to those of "normal" spi- 
ral galaxies of the same stellar mass. Figure 2 of Bigiel 
et al. (2008) shows that in spirals where H2 dominates 
in the central regions, Esfr decreases from 0.01-0.1 Mq 
yr -1 kpc -2 near the galaxy centers to less than 0.001 
Mq yr _1 kpc -2 near the edges of the optical disk. The 
SFR surface density profile of GASS35981 is much more 
similar to those of Hi-dominated dwarf irregular galax- 
ies (Figure 3 of Bigiel et al.), but these galaxies have 
considerably smaller stellar masses. 

In the right panel of Figure HI we plot EW(Ha) as a 
function of position along the slit. Since the equivalent 
width is the ratio of emission line flux to the underlying 
continuum, it is an excellent proxy for ratio of current 
star formation to past-averaged star formation (i.e., pre- 
existing stellar mass), conventionally referred to as the 
6-parameter. The exact relation between the two de- 



http://www.mpa-garching.mpg.de/SDSS/ 



Another possible worry is that b parameters and star 
formation rates measured through small apertures are 
not representative of the galaxy as a whole; however, 
we obtain nearly identical trends with radius when we 
estimate b and Ssfr from photometry that has been 
azimuthally-averaged in elliptical annuli across the whole 
galaxy. 

To test whether the age of the stellar population in 
the outer disk is consistent with the short formation 
timescale deduced from the Ha equivalent widths, we 
examine the radial trends in stellar population age as 
measured by the D4000 n index (Figure [7]). Although 
£>4000 n is known to vary with metallicity as well as age 
(Kauffmann et al. 2003), the metallicity of GASS35981 
appears to be fairly uniform within the restricted range 
of radii over which we measure it (with only one point in 
Figure [7] at significantly sub-solar metallicity; see § 4.2), 
and so for this radial regime D4000„ is the highest-S/N 
proxy for stellar population age that we can measure. 

We see in Figure [7] that Z?4000„ appears to decrease 
from ~ 1.7 in the center of the galaxy to 1.3 — 1.4 beyond 
~ 10 kpc. For a simple single stellar population (SSP) 
model from Bruzual & Chariot (2003), L>4000„ = 1.3 
indicates a stellar population age of 500 Myr, while 
-D4000 n = 1.7 occurs at an age of ~ 3 Gyr. These charac- 
teristic timescales grow larger for more realistic extended 
star formation histories; this will be discussed further in 
§ 5. The S/N of our spectrum is not high enough to mea- 
sure -D4000 n at very large radii, but the fact that such 
young ages are implied for R = 10 kpc is suggestive that 
the outer disk may be even younger. 

4.1.1. Analysis of the CO Measurements 

Given the order of magnitude difference between the 
measured HI mass and the measured molecular gas mass, 
it is important to determine whether the latter is consis- 
tent with the derived star formation rate. Bigiel et al. 
(2008) show that SFR surface density is directly pro- 
portional to H2 mass surface density according to the 
relation: Esfr = 10~ 31 Eh 2 (in units of M© yr _1 kpc -2 
and Mq pc -2 , respectively). This relation holds over a 
large range in Eh 2 (3-50 M Q yr _1 pc -2 ) and does not 
appear to vary as a function of any larger scale property 
of the galaxies in the sample. 

When comparing SFR and H 2 mass for GASS35981, 
we must carefully account for possible aperture effects, 
because the SFR surface density is measured through a 
narrow slit, while the H2 masses are estimated in three 
much larger apertures with FWHM of 22". However, we 
have seen that in GASS35981 the SFR surface density ap- 
pears to be quite uniform across the whole galactic disk. 
The CO line flux measured in our central pointing corre- 
sponds to an average H2 surface density across the beam 
of Eh 2 = 3.6 ± 0.1M Q pc -2 , under our assumed conver- 
sion factor Xco = 4.4M /L'. Adopting the Bigiel et al. 
relation between Eh 2 and Esfr yields a predicted SFR 
surface density E SFR = 0.0029 ± 0.0001M Q yr" 1 kpc" 2 , 
which is in excellent agreement with the observed Esfr 
in this galaxy. 

In the two offset pointings where we do not detect 
any CO, our most stringent upper limit (in the North- 
ern pointing) implies that Eh 2 < 2.4M pc -2 , and hence 
Esfr < 0.0019. This is 50% lower than the median Esfr 



we measure from our spectra, though from Figure |H] we 
see that our photometric estimates of Esfr exhibit a dip 
in the region covered by the northern pointing, which 
could explain the low H2 mass. 

Alternatively, variation in the conversion factor Xco 
could explain the discrepancy; since below we will present 
evidence that the metallicity is low in the outer regions 
of GASS35981 (which implies a lower amount of CO per 
unit of H 2 ), this possibility seems plausible. We note 
that Bigiel et al. (2008) neglected variations in Xco m 
deriving their SFR-H2 relation, and so the ~ 0.3 dex 
scatter they measure likely includes a contribution from 
variations in metallicity /Xco- Though the single dis- 
crepant point here is too little to draw conclusions from, 
in a future paper we will examine possible variations in 
Xco across the COLD GASS sample. 

4.1.2. Total Star Formation Rate 

Finally, we estimate a total star formation rate for 
GASS35981 by extrapolating the SFR surface density 
measured for each spectral bin along the slit to the an- 
nulus that spans the same radial range as that portion 
of the slit, and that extends half-way around the galaxy, 
where it meets the corresponding annulus extrapolated 
from the other side of our slit. The aperture corrections 
are large, but it is instructive to compare the resulting to- 
tal, dust-corrected SFR with our purely photometric es- 
timate. We find that SFR Hq = 3.7±O.3M yr _1 , which 
is slightly higher than the value of 3.O±O.4M yr _1 esti- 
mated from the integrated photometry using the method 
described in Schiminovich et al. (2010). We find no evi- 
dence for any additional obscured component of star for- 
mation in GASS35981, as there is no counterpart in the 
Faint Source Catalog of IRAS (Moshir et al. 1992), to a 
limit of 0.2Jy at 60/im, which corresponds to a luminosity 
ofi 6 o~lO fo L . 

4.2. Metallicity Gradient 

We utilize the [Nil] 6584 to Ha emission line ra- 
tio as our primary metallicity indicator. We estimate 
metallicity from the relation of Pcttini & Pagel (2004): 
12 + \og(0/H) = 9.37 + 2.03 x N2 + 1.26 x A2 2 + 
0.32 x A2 3 where N2 = log([NII]6584/Ha). The re- 
sulting mctallicities are plotted as a function of position 
across GASS35981 in Figure [8] In the bulge and inner 
disk of the galaxy (R < 15 kpc), GASS35981 exhibits 
near-solar metallicity. However, at large radii, we see a 
sharp drop in metallicity of ~ 0.5 dex, at a radius corre- 
sponding to the start of the blue, faint outer disk of the 
galaxy (Figure [3]). This is also the region of the galaxy 
where the b parameter exceeds unity. 

Sharp breaks to low metallicities are not common in 
the portion of the larger GASS sample with follow-up 
MMT spectroscopy that has been analyzed so far. We 
find an average gradient in \og{0/H) of —0.14 dex/r e , 
compared to —0.41 dex/r e for GASS35981. Less than 
10% of the galaxies in our sample exhibit metallicity be- 
low 12+log(0/H)< 8.5, and in most of these cases, the 
low metallicities are seen over the whole galaxy. Only 
~ 15% of the 62 galaxies analyzed so far have metallic- 
ity gradients as strong as GASS35981, and even fewer of 
these exhibit a sharp drop. We note that there is one 
well-known local galaxy with a metallicity gradient sim- 
ilar to GASS35981: M101 (van Zee et al. 1998). The 
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Figure 8. Top: Metallicity of GASS35981 as a function of posi- 
tion along the slit, estimated from the [NII]/Ha emission line ratio 
(Pettini & Pagel 2004). Solar metallicity is marked with the dotted 
line. Bottom: Similar metallicity estimates including the ratio of 
[OIIIJ/H/3, also from Pettini & Pagel, as described in the text. 

importance of this similarity is unclear, however, since 
the two galaxies differ in other respects (e.g., the distri- 
bution of star formation; Kuntz & Snowden 2010). 

As a check on our [Nil] /Ha metallicity esti- 
mates, we also compute metallicity using 03N2 = 
log([OIII] 5007/H/3/QNII] 6584/Ha)) (Pettini & Pagel 
2004) , which is plotted in the bottom panel of Figure [SJ 
Though intrinsically more accurate, this indicator pro- 
vides somewhat poorer constraints for our galaxy, be- 
cause H/3 is not always detected at > 3er in our spectra. 
We note, however, that metallicities measured in this way 
are entirely consistent with those from [Nil] /Ha , as can 
be seen by comparing the two panels in Figure[5] We also 
use ([SII] 6717+ [SII] 6731)/Ha as an additional check, 
as this line depends on metallicity in much the same way 
as [Nil] (Dopita et al. 2006). Although it is detected 
with lower S/N in our spectrum, and is more affected by 
sky subtraction residuals due to its redder wavelength, 
the [SII] /Ha ratio exhibits a similar drop at large radii. 
This indicates that peculiar Nitrogen abundances are not 
responsible for the effect that we see. 

Finally, in Figure H] we plot metallicity as a function of 
EW(Ha). There is clearly a strong correlation: regions 
of the galaxy with the highest ratio of current to past- 
averaged star formation rate are also the most metal- 
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Figure 9. Metallicity from the [NII]/H« ratio as in Figure[8] but 
plotted as a function of EW(Ha). Increasing Ha EW is strongly 
correlated with decreasing metallicity in GASS35981. Solid line in- 
dicates the relation followed by normal star-forming galaxies drawn 
from SDSS, as described in the text. 

poor. We also over-plot the best-fit relation between 
EW(Ha) and [Nil] -based metallicity for star- forming 
galaxies drawn from the SDSS MPA/JHU value-added 
catalog. These measurements apply to the region of the 
galaxies sampled by the 3"-diameter SDSS fiber. The re- 
lation between star formation rate and metallicity mea- 
sured within GASS35981 follows that exhibited by the 
general population of star-forming galaxies in SDSS. We 
note, however, that the galaxies with high EW(Ha) and 
low metallicities are generally much lower mass systems. 

5. DISCUSSION 

We have seen that GASS35981 is currently undergoing 
star formation that is evenly spread across the galaxy. In 
the outskirts, however, a sharp drop in metallicity and 
a correspondingly large increase in the EW(Ha) suggest 
that the bulk of the stellar mass in this region has been 
formed only recently, and out of a low metallicity reser- 
voir of gas. 

In this section, we ask how long ago the gas began 
forming its stars, and we consider the processes that may 
have delivered the gas to the galaxy. 

5.1. A Model Star-formation History 

We construct a simple toy model of the galaxy's star 
formation history (SFH), in the bulge (R < 5 kpc), inner 
disk (5 < R < 15 kpc) and outer disk (R > 15 kpc) of the 
galaxy. Using constraints derived from the galaxy's cur- 
rent star formation rate and stellar mass surface density, 
we build a simple model that is consistent with observa- 
tions in all three locations. 

Our observations show that the current star-formation 
is uniformly distributed across the galaxy, so we postu- 
late that the overall SFH can be modeled by a period of 
constant star formation of some length that is the same 
everywhere, superimposed on an older stellar population 
that dominates at the center and is virtually absent at 
large radii. 

Our model incorporates the following constraints: 

• At R > 15 kpc, the metallicity is low and the 
b parameter is extremely high. Let us assume 
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Figure 10. Schematic of a model star formation history that 
replicates key features of GASS35981's stellar population in each 
of three radial regimes. We model the current star formation as a 
single episode of length 1 Gyr proceeding at the same rate every- 
where across the galaxy (normalized to SFR = 1 in this diagram). 
The underlying old stellar population is modeled as another 1 Gyr- 
long burst occurring at z = 1, with amplitude constrained by the 
observed present-day stellar mass surface density in each radial 
regime. 

that this indicates all of the existing stellar mass 
has been built during the current star-formation 
episode. For the three outer spectral bins with the 
lowest metallicities, the timescales for forming the 
entire stellar population at the currently observed 
star formation rate range from 0.7-2.0 Gyr. We 
thus adopt 1 Gyr as our reference value. 

• We then subtract the stellar mass formed in the 
current episode from the observed stellar mass in 
the bulge and the inner disk, and assign the re- 
sulting mass to the old stellar population. In the 
inner disk, the mass surface density in the old stel- 
lar component is ~ 10 times that in the outer disk, 
while in the bulge region, this factor is closer to 
60. For simplicity, we will concentrate all of the 
star formation needed to build this much mass into 
a single episode of 1 Gyr length ending at z = 1 
(7 Gyr in the past). 



A simple schematic of this model SFH is shown in Fig- 
ure [TUJ We then use Bruzual & Chariot (2003) models 
to predict Z34000 n in the inner disk and bulge, obtaining 
values of 1.3 and 1.5, respectively. These values are in 
good agreement with the measured values of 1.4 and 1.6 
in the corresponding radial bins. Thus, through only a 
simple partition of mass into old and young components, 
we can reasonably reproduce one of the key spectral fea- 
tures of GASS35981. We note that varying the formation 
time of the old component or doubling its timescale to 
2 Gyr changes the results very little. Likewise, adjust- 
ing the length of the more recent star-forming episode 
within the range 0.7-2 Gyr does not significantly affect 
our results. 

Finally, we can estimate the total fraction of the stellar 
mass added to the galaxy in the most recent star forma- 



tion episode and find that it is around 20%. Note, how- 
ever, that at its current star formation rate, GASS35981 
will not exhaust its HI reservoir for another 5-7 Gyr. 

5.2. Was the Gas Captured from Another Galaxy? 

We have concluded that a new episode of star forma- 
tion began approximately 1 Gyr ago in GASS35981. It 
is tempting to speculate that renewed star formation in 
this galaxy is connected to the acquisition of its large HI 
reservoir. 

So how did GASS35981 acquire its gas? One possibility 
is that the gas was accreted as the result of an interaction 
with another gas-rich system. GASS35981 does have two 
luminous neighbors at projected distances of 300 kpc and 
400 kpc with redshift differences of < 400km s _1 — i.e., 
GASS35981 does appear to be a member of a loose group. 
We do not, however, see any signs of an interaction or 
merger in either the rotation curve or optical morphology 
ofGASS35981. 

Furthermore, a scenario in which GASS35981 is left 
with 2 x 10 10 M Q of HI after an encounter with either of 
its neighbors seems fairly implausible. Simulations pre- 
dict that no more than 20% of the gas mass of any donor 
galaxy should be stripped in an encounter (Bournaud, 
private communication), so we would require a donor 
galaxy with > 10 11 M Q of HI! Such an encounter could, 
however, be responsible for jostling an already-present 
HI reservoir out of equilibrium, causing it to form stars. 

If the gas was not acquired from another passing 
galaxy, one might speculate that GASS35981 accreted 
its HI reservoir directly from the surrounding intergalac- 
tic medium at some point in the past. Since blind HI 
surveys such as ALFALFA (Giovanelli et al. 2005) and 
HIPASS (Barnes et al. 2001) do not find dark HI clouds 
of 1O 1O M , the gas must have entered GASS35981 from 
an ionized phase if it entered all at once or over a short 
timescale. Since its two identified neighbors suggest that 
GASS35981 resides in a group, the presence of intra- 
group gas could be fueling an unusually high accretion 
rate onto GASS35981. Indeed, accretion from the intra- 
group medium has been speculated to be responsible for 
a number of other peculiar star- forming systems (e.g., 
Beaulieu et al. 2010). 

Alternatively, the gas reservoir could have been built 
slowly through multiple accretions of smaller gas clouds 
or streams, which could be either neutral or ionized. 
Under this scenario, star formation would need to be 
suppressed somehow during the buildup of the reservoir. 
The galaxy formation models of Birnboim et al. (2007) 
exhibit quiescent, reservoir building periods similar to 
what would be needed here, but in general they apply 
to somewhat higher mass galaxies than GASS35981, and 
also may not be valid at z ~ 0. Multiple minor merg- 
ers with gas-rich dwarfs could also supply the gas, but 
it becomes even harder in this case to imagine how the 
HI could build up over time rather than form stars with 
each new accretion event. 

6. CONCLUSIONS 

We have reported on the remarkable galaxy 

GASS35981, a disk galaxy with stellar mass 

M* = 2 x 10 10 M which contains an additional 

2.1 x 10 10 Mq of HI gas. Millimeter observations 



11 



indicate a molecular gas mass only a tenth this high. 
Through follow-up long-slit spectroscopy, plus SED 
fitting using our UV through optical photometry we 
have shown that: 

• Star formation is evenly spread across the 
galaxy, at a surface density of Ssfr = 
0.003 M yr- 1 kpc~ 2 . 

• The proportion of the stellar mass contributed by 
the current star-formation episode rises towards the 
outer regions of the galaxy, reaching a peak at R = 
30 kpc, where the entire disk must have formed in 
the past Gyr. 

• Interstellar metallicities exhibit a sharp drop at 
R > 15 kpc, coincident with a sharp rise in 
EW(Ha). This is consistent with recent infall of 
lower- met allicity gas (Tinsley & Larson 1978). 

• The Ha rotation curve is regular and symmet- 
ric, and reveals no signs of a recent interaction or 
merger that could have deposited the gas and/or 
triggered the recent star formation episode. 

The main conclusion from our observations is that 
GASS35981 appears to be in the early stages of forma- 
tion of its outer stellar disk. 

We are not able to provide conclusive answers to ques- 
tions pertaining to the origin and fate of the gas in 
this galaxy with this data set alone. Scenarios in which 
the gas was acquired in a recent merging event are dis- 
favoured because of the extremely regular kinematics of 
the disk. The HI mass of GASS35981 is too large to be 
easily explained by gas transfer from a passing galaxy. 
We therefore speculate that GASS35981 acquired its gas 
directly from the intergalactic medium. Although our 
observations show that the stars in the outer disk formed 
within the last Gyr, this does not mean that the gas was 
also acquired less than 1 Gyr ago. It is also unclear 
whether GASS35981 will continue forming stars in its 
current low-efficiency state, or whether the gas will flow 
inwards towards the bulge, and GASS35981 will even- 
tually develop into a more normal massive spiral galaxy 
with a star formation surface density that decreases as a 
function of radius. 

Questions concerning the eventual fate of the gas can 
be addressed using the larger samples that will be pro- 
vided by the full GASS and COLD GASS surveys in the 
future. By studying trends in SFR surface density, mean 
stellar age, metallicity, and stellar mass profiles as a func- 
tion of atomic and molecular gas content for complete 
samples of galaxies, we hope to map out evolutionary 
sequences in disk galaxy formation. 

Answers to questions concerning the origin of the gas 
will likely require a different approach. Our compari- 
son of the HI linewidth of GASS35981 with its CO line 
width and Ha rotation curve yield tantalizing hints that 
the atomic gas may not be in equilibrium with the rest 
of the galaxy. In addition, the HI spectrum in Figure [1] 
is clearly asymmetric about the line center. High reso- 
lution HI mapping of GASS35981 will be needed to un- 
derstand the dynamical state of the gas in more detail. 
Even so, such observations are unlikely to prove that the 
HI originated from a more diffuse (and unseen) reservoir 



of IGM gas. This can only be done if we are able to 
find tracers of this gas, for example absorption lines in 
the spectra of background quasars that arise when the 
quasar light passes through the circumgalactic medium 
of the galaxy (Cen & Ostriker 1999). These tracers must 
then be linked with galaxies like GASS35981. 
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